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Abstract 

The general formula for the interaction potential between two point electric charges is derived. 
This analytical formula has the correct asymtotic behaviour at large distances between two inter- 
acting charges. The derivation of this formula is based on the closed analytical expression for the 
p I Uehhng potential obtained earlier (A.M. Frolov and D.M. Wardlaw, ArXiV: 1110.3433 [nucl-th]). 
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In classical electrodynamics the interaction V between two point electric charges qie and 

2 

q2e is described by the Coulomb law, i.e. V ~ ^^f^- This expression for V does not change 
its form in the non-relativistic quantum mechanics [lf| , when the Planck constant h is finite. 
However, if we also assume that the speed of light c is finite, then Quantum Electrodynamics 
leads us to the necessaty to modify the Coulomb law. The main correction is related to 
the vacuum polarization in the spatial areas close to the electric charges. In the lowest 
order approximation such a correction is represented by the Uehling potential which is 
correct at short interparticle distances. However, the Uehling potential does not provide 
the correct long-range asymptotics, i.e. its behaviour at large r is not correct. Therefore, it 
must be modify to represent the interaction between the two point electric charges correctly 
at arbitrary, in principle, interparticle distances. This problem has substantial scientific 
and methodological interest. In particular, the Coulomb potential modified by a number 
additional terms, which represent small QED-corrections, can be used directly in highly 
accurate calculations of bound states in few-electron atoms and ions. On the other hand, 
it is interesting to see the transformation of the usual Coulomb potential under the rules of 
modern Quantum Electrodynamics. 

In our earlier study [3] we have derived the closed analytical expression for the Uehling 
potential U{r) which describes the lowest order QED correction (upon the fine structure 
constant a) to the regular Coulomb potential. This correction is related with the polarization 
of the vacuum produced by a point electric charge. In jsl] it was shown that the total 
interaction potential is written in the form (in atomic units h = 1, nie = 1 and e = 1) 
, Qr^ 2a r/^ \ /2r\ r /2r\ /5 \ /2r 
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where f/(r) is the Uehling potential |2| and a = |^ = c^^ is the fine structure constant. Here 
and everywhere below in this study the notation h stands for the reduced Planck constant 
[fi = ^), e designates the absolute value of the electric charge of the electron and c is the 
speed of light in vacuum. In Eq.([T]) the function KQ{a) is the modified Bessel function of 



zero order (see Eqs. (8.432) and (8.447) from 
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Ko{z) = exp(-zcosht)rft = ^(V'(A; + 1) + ln2 - Inz) , 
Jo ^0 2^''{k\y 

where ip{k) is the Euler psz-function defined by Eq.(8.362) from ^. The functions Kii[z) 
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and Ki2{z) in Eq.([T]) are the recursive integrals of the Kq{z) function, i.e. 

roo poo 

Kii{z) = / Kio{z)dz , and Kin{z) = / Kin-i{z)dz , (2) 
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where n > 1 and KiQ{z) = Kq{z). 

The formula, Eq.([l]), is corret mathematically, but it has a wrong asymptotic at large 
distances r, e.g., at r > lOaao, where ao = ^^^-^ is the Bohr radius. First, note that the 
short range asymptotic of the (j){r) potential, Eq.(Il]), takes the form 

Q (-1 : r ^ 
~3 

where 7 ^ 0.5772156649 ... is the Euler constant (see, e.g., |4|]) and Ina ^ -4.92024365857. 
The long-range asymptotics of the potential (f){r), Eq.([T]), is (in atomic units) 
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<P{r) |,^^^^{l + ^l^exp(--r)} (4) 

This means that the long-range asymptotics of 0(r) decreases with r exponentially. This 
has no physical sense, since the exponential function vanishes at r — )■ +00 very rapidly. The 
corresponding correction to the potential at large distances will always be overweighted by 
the lowest order QED correction to the electromagnetic field (or EM- field, for short) (E, H), 
or to the pure electric field E in our case, when H = 0. This correction provides power-type 
dependence at large distances and it is directly related to the non-linearity of the Maxwell 
equations for the EM-field j^l, jst. 

Let us obtain the lowest order QED correction to the electric field E by considering one 
point electric charge e which is at rest in the center of coordinates x = 0, ?y = and z = 0. 
The corresponding Lagrangian L which generates the lowest order QED correction to the 
energy of the electric field is written in the form (in regular units) (see, e.g., (sl) 

L = 1e^ + ^J^E^ = -(V0)' + ^J^(V(^)^ = LU>{r)\ (5) 

where E = — V(/)(r). By varying the potential 0(r) in the following equation (the funda- 
mental equation of least action p) 

5 J L[(j){r)ydr = (6) 

one finds the differential equation 

d(f) e^h /d(j)\3 C , , 



where C = — Assuming that (/)(r) = ^ + ip{r), where '0(r) is a very small correction 
we reduce the last equation to the form 



e 



From this equation one finds 



Therefore, the total interaction potential 0(r) is 
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In the relativistic units h = 1, c = 1 and a = ^ the last formula takes the form 
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This expression is the sum of the Coulomb potential and small correction which exactly 
coincides with the Wichmann-KroU potential jsj. In atomic units the expression, Eq.( !TT|) . 
takes the form 



r = - 
r 



Qr 2g^a^ Q 2Q^a^ 



2257rr4J r 2257rr5 
and Wichmann-KroU potential Wk[t) is 



(12) 



Note that this potential is not regular at the origin. However, as follows from the Appendix 
such a singularity is formal, since Eq.(I7]) is not correct at very small and very large r (see 
Appendix). The singularity at small r can be removed by using substitution r — )■ r + a. Such 
a substitution removes singularity of the VFft-(r) potential at r = and changes its behaviour 
at short interparticle distances r < a. However, at such distances the Wichmann-KroU 
correction is not important and its contribution is significantly smaller than contribution 
from the Uehling potential. On the other hand, at large distances, i.e. for r > lOaao (and 
even for r > aag) the contribution from the Wichmann-KroU potential is substantial and 
cannot be neglected (the correction from Uehling potential vanishes exponentially with r). 
Finally, we can write the correct expression for the Wichmann-KroU potential Wk{t) 

Mr) = WK(r) = 7^ " , (14) 
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Note that this potential is always negative. In general, at r ^ aQ this potential is 
very small by its absolute value. However, its overall contribution rapidly increases with 
the nuclear charge Q. Furthermore, it decreases with the distance r as ~ r~^, i.e. non- 
exponentially. It is clear that at large distances the Wichmann-Kroll correction will always 
overweight the contribution form the Uehling potential U{r), Eq.([T]). Thus, the following 
interaction potential (in atomic units) of the two point electric charges 

= « + U(r) + H/,M = « + . [(1 + ^)A'„(?^) - f 



3a^^ ^ a ^ 6a; ^ a 
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has the correct asymptotic behaviour both at small and large interparticle distances r. In 
the case of interaction between two point electric charges gie and q2e we need to replace in 
Eq.f lT5|) the factor Q by the product qiq2. The potential $(ri2) takes the form (in atomic 
units) 

$(^^2) = ^ + ^^^^^^ ■ [fl + Tk.\Ko(^) - —Kiii^ 
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3a^^ ^ a ^ 6a ^ a 
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where ri2 is the distance between particles 1 and 2. In regular units one needs to replace 
in Eqs. ffTSj) - (fT6|) the dimensionless value a of the fine structure constant by the factor 
a-^. This produces the regular expression for the Wichmann-Kroll potential for atomic 



systems. Regularization of this term for systems of heavy particles requires an additional 



investigation. 

Thus, we have derived the formula for the interaction between two point electric charges 
which includes the lowest order QED corrections to the Coulomb potential. Our formulas, 
Eqs.f llSp - f|T6|) . provide the correct asymptotic behaviour at arbitrary interparticle distances. 
The formulas, Eqs. (fT5|) - f|T6|) . can directly be used in highly accurate computations of the 
bound states in few-electron atoms and ions. The explicit computation of the electric field 
strength E, i.e. the spatial gradient from the $(ri2) potential, Eq.(IT6l). is straightforward, 
but the Fourier resolution of the corresponding 'electrostatic' field is extremely complicated. 
Note also that for the potential $(r), Eq. ffTSj) . the condition A$(r) = AirQdlr), where 
A = V ■ V = div{grad...), is not obeyed. 
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Appendix 



The equation, Eq.([7]), is a qubic equation in respect to the radial derivative of the unknown 
potential (l){r). It can be re- written in the form 



y +y + 
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where y = As follows from the general theory of qubic functions this equation has only 
one real root (the discriminant of this equation is negative). In fact, this equation is easily 
reduced to the form of a monic trinomial (see, e.g., [TJ) 



(18) 



where p 



and q 



The Cardano method gives the only real root of Eq. ([18 
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In other words, we obtain the following differential equation 



dr 
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At r — i- oo this equation has the solution 0(r) = 2(|)2 ■ r + c, where c is some numerical 
constant. At r — > the analogous solution takes the form 0(r) = — 3^^ ■ + c, where c is 
some numerical constant. Both these asymptotics do not have any physical sense. Therefore, 
our procedure used in the main text to approximate the exact solution of Eq.([7]) is correct 
only for some restricted interval of r— variation. In other words, the substitution r — )■ r + o; 
used in the main text to regularize the Wichmann-Kroll potential at small r is correct. 
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